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Abstract 

A new kind of tridimensional scalar optical beams is introduced. These 
beams are called Lorentz beams because the form of their transverse pat- 
tern in the source plane is the product of two independent Lorentz func- 
tions. Closed-form expression of free-space propagation under paraxial 
limit is derived and pseudo non-diffracting features pointed out. More- 
over, as the slowly varying part of these fields fulfils the scalar paraxial 
wave equation, it follows that there exist also Lorentz-Gauss beams, i.e. 
beams obtained by multipying the original Lorentz beam to a Gaussian 
apodization function. Although the existence of Lorentz-Gauss beams can 
be shown by using two different and independent ways obtained recently 
from Kiselev [Opt. Spectr. 96, 4 (2004)] and Gutierrez- Vega et al. [JOSA 
A 22, 289-298, (2005)], here we have followed a third different approach, 
which makes use of Lie's group theory, and which possesses the merit to 
put into evidence the symmetries present in paraxial Optics. 

1 Introduction 

Optical beams are electromagnetic fields with a well distinguishable mean direc- 
tion of propagation (that, from now on, we indicate as z axis), in the nearness 
of which the most part of field's energy is contained during propagation. After 
the publication of the fundamental work due to Durnin et al. J^, issued in 1987 
about non-diffracting beams in free-space, there was a certain number of sci- 
entifical investigations on optical fields with the structure of beams and which 
possess, or approximate at least, the diffraction features of the aforesaid beams 
Indeed, as an ideal non-diffracting beam is physically unrealizable, be- 
cause on each plane z = constant it carries an infinitive amount of energy, it 
is possible to obtain diffraction-free like fields only with approximation: these 
fields are also known as pseudo non-diffracting beams. Examples of that are 
the well-known Gaussian beams, Bessel-Gauss beams 4 or other optical fields 
recently introduced by casting the propagation problem in coordinate systems 
different of rectangular and circular ones 5 that share the properties of man- 
taining a non-diffracting behaviour only inside a limited spatial range named 
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Rayleigh distance. In the present work, we wish to introduce another class of 
pseudo non-diffracting rcahzable beams that we wiU caU Lorentz beams {LB for 
short) as well as their Gaussian apodizated version, that, to authors's knowledge, 
were never been study before today. If, from a theoretical point of view, the re- 
search of new kinds of optical beams is interesting, it is particularly stimulating 
in the present case because of the physical realizability of the proposed field. 
This realizability is not due to a Gaussian term, as usually happen for almost all 
other known optical beams, but from a practical point of view it was shown |7]- 
[Hj that certain laser sources produce fields that shows fundamental variations 
with respect to the canonical Gaussian beam. As well-known Gaussian beam is 
a minimun uncertainty field i.e. it possesses the minimum achievable angular 
spreading once the spatial extension is fixed; for certain laser sources, e.g. dou- 
ble heterojunction (DH) Gai-^xAlxAs lasers, which produce highly diverging 
fields, a Gaussian description for the transverse fields fails: in this case it was 
shown that a Lorentzian distribution is a better approximation, as it take into 
account of the higher angular spreading, being equal the spatial extension [Hj. 
The paper is organized as follows: in next section we introduce the LB and we 
study their propagation under Fresnel or paraxial approximation. In particular 
we give the closed-form expression for these kind of fields on a generical (x, y) 
plane. Furthermore, we will make use of a theory group approach to introduce 
another class of optical beams obtained by multlipying a Lorentz beams with a 
two-dimensional Gaussian envelope. 



2 Lorentz Beams and paraxial propagation 

Let us suppose we have on a source plane, that we decide to be {x,y,z = 0) 
plane, the following scalar field distribution 

A 1 1 

^"("^'^^ ^ WxWy [1 + [x/wxy] [1 + iy/wyr] 

where A is a constant value and Wx and Wy are parameters related to the beam 
width, with A, Wx and Wy € 3?. This kind of field is the product of two functions 
of X and y variables which have the form of a Lorentzian function of parameter 
Wx and Wy. The Lorentzian is a well known bell-shaped curve used principally 
to describe the spectral lines of simple dinamical systems usually present in 
physics. Starting from the field in we wish to derive the form that such 
a field acquires during free propagation. To do this, we pass to the Fourier 
domain, calculating the plane waves spectrum on z = 0. We have 

f f A 1 1 

WxWy [1 -I- (X/Wx)^] [1 -I- {y/Wy)'^] 

X eKp(—2n:qy)dxdy] 
A f 1 

= / r-i , / / — exp{~-2inpx)dx x (3) 

WxWy [1 + {x/Wx)^\ 

exp{—2inqy)dy 



[1 + (y/wy) 

= An'^ exp{~2TT\p\wx) exp{—2Tr\q\wy) (4) 
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where AQ{p,q) is the complex spectrum, p and q are the conjugated variables 
to X and y variables respectively. Once we know the spectrum of plane waves 
on z = we can calculate it on a generical (x, 2/)-plane as follows 

Az{p, q) = Ao{p, q) exp(2i7rmz) = Att^ eyij)[—2T:\p\wx) exp(— 27r|q|z«y) exp(2i7rmz) 

(5) 

because it is well known in which fashion a plane wave propagates in free-space. 
The parameter m is a complex variable that must fulfil the relationship 

m = v/(l/A2-p2_g2) (6) 

and it is imaginary for evanescent waves, in which p^ -\- q^ > 1/ , and real for 
homogeneous waves where p^ + q^ < l/A^, where A is the wavelength. We know 
that a field is said to be homogeneous when the spectrum A^ (p, q) is different 
from zero only inside the circle p^ + q'^ < l/A^. 



2.1 Fresnel or paraxial limit 

Observing the Q we deduce that if the values of Wx and Wy are enoughly 
greater than wavelength A, the majority contribution to the field arises from 
homogeneous waves having the amplitude Az(j),q) corrisponding to points, in 
Fourier plane (p, q), near to the origin. In this case one is authorized to do the 
following expansion 



,,1 2 2^ 1 (p^ + g^)A 
n.= ]J{^-p'-q')^j 2 

by retaining only the first two terms in the series, so that the plane wave spec- 
trum in ^ becomes 

^parax^p^ g) = Att'^ exp{~2Tr\p\wx - inXzp'^) exp(-27r|g|wy - inXzq'^) exp{ikz) 

(8) 

in which k = 2tt/X is the wave number. Now, as we know the plane wave 
spectrum on z we can also derive the full form of the field if we are able to 
inverte the two-dimensional Fourier transform 

Vix,y,z) . / /[A.^exp(.Mexp(-2.b|..-..A.p^)x 

X ex.p(^27T\q\wy — i7TXzq^)]dpdq 

^ An^eM^kz) f eM-2Mp\^.-^^X^P^)eM^2np)dp. 

J oo 

X / ex-p{—2n\q\wy ~ iirXzq^) exp{i2TTq)dq (9) 

J oo 

It is important to note that the role of paraxial approximation was to give a plane 
wave spectrum A^"'""^ (p, g) factorized in two terms, each one depending only 
from a single Fourier variable, p or q; this is a properties that was not fulfilled 
by the exact spectrum in 101 for the presence of the term exp(i27rmz). As a 
consequence of that also the complex field V{x,y,z) is in a similar factorized 
form. To obtain a solution of jnj let us focus on the integral 

I — cxp{—2'k\p\wx ~ iirXzp'^) exp{i2Trp)dp (10) 

J oo 
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(11) 



It can be written as 

/"OO 

I = exp{—2npwx — iTrXzp^) exp{i2'Kpx)dp + 

Jo 

nO 

+ / exp{2TrpWx — iirXzp^) exp{i2TTpx)dp (12) 
For the first of two integral in right-hand side we have 

/•OC 

I = cxp{—2'KpWx — iirXzp^) exp{i2TTpx)dp = 

Jo 

= exp(c) / exp[~{ap + b) ]dp ^ — / exp(~s^)ds 

Jo a Jb 

= ^^^[ r exp{-s^)ds - [ exp{-s^)ds] 
a Jo Jo 

= ^g ^xp[-K -^)V»Az] _ ^^^^^^^^ _ (^3) 

where we have defined the foUowing auxiharies variables 

= inXz (14) 

b = Tr{wx — ix) / a (15) 

c = b^ (16) 

s = ap + b (17) 

and erf{x) is the usual error function erf{x) = (21^) f^^ exp{—s^)ds 6 . The 
other integral appearing in (|ll|l can be easily calculated by observing that it is 
equal to that just derived in after having substituted the variable x with 
—X. On utilizing this result we finally obtain the full form of the field 

Vix,y,z) = ^^^^^^[v+{x,z) + Vx-ix,z)][V^{y,z) + V-{y,z)ll8) 
where 

4_ , , exT>\Tr(wr ± ir)'^ /iXz] , , , , i — ^ — ,^ 

V^(r,z) = ^ J-^ 41 - er/[7r(u;^ ± ir)/\/i7rAz]} (19) 

and r — x,y. Equation IjlSjl is the principal result of the present work and in 
next sections we analyse more in detail the propagation features of these beams. 



2.2 Propagation and diffraction- free range 

We expect that the field in H18() changes its shape during propagation as a 
consequence of diffraction. It is well known, however, that it is possible to 
define a diffraction-free range, i.e. a linear distance on z-axis, under which the 
beam remains essentially unchanged. To do this we write the beam to values 
near to the source plane {z = 0) by taking advantage of the expansion of error 
function for high values of its argument 

erf{s) « 1 + ^ exp(-s^) g(-l)'=ItiM (20) 

77 ^ — ' S 

k=l 
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with \s\ >> 1, — 7r/2 < arg(s) < tt/2 and T is the gamma special function 
On utilizing this expansion it is easy to verify that (|18() reduces to Indeed 
one obtains, by keeping only the first term 

Vix,y,0) = ^lim-^[V+{x,z) + V-{x,z)][V+iy,z) + V-{y,z)] = 

4 z^O l\z " " 

= — — hm — ^ , + ^ ' ' X 

4 2^0 lAz y/inr[Wx ~ ix) y'Tnr(Wx + ix) 

r(l/2)VlAi ^ ^(l/2)^^Ai j_ 
V7r7r(wj^ - iy) ^/TTTr{wy + iy) 
A-K^ , 1 1 1 1 , 



-][- 



4 Wx ~ ix Wx + ix Wy — iy Wy + iy 
A 1 1 

WxWy [1 + (X/Wx)'^] [1 + (y/Wy)'^] 



(21) 



where we used the identity r(l/2) = y/n. The shape of the field will be practi- 
cally unchanged as far as the second term in expansion will be negligible respect 
to the first. This conducts us to the following condition, 

^(1/^) » m^) (22) 



\Tr{wr ± «r)/A/i7rAzp \n{wr iir)/\/inXz\'^ 
where r = x,y. Equation leads to 

r(l/2) >> r(3/2)^z (23) 

and finally (recall that r(3/2) = v^/^) we arrive to the evaluation of the 
diffraction-free range (or Raileigh distance) for this kind of beams 

ZR = (24) 

In H23|) we have let Wx ^ Wy ~ w to simplify the analysis. If one does not 
make this assumption there will exist two different diffraction-free ranges, one 
to a;-axis and another to y-axis. In figure we report the modulus of field's 
amplitude evaluated to different distances from the source. In particular we 
observe that, for distances sufficiently shorter than the effect of diffraction 
is neagligible, as expected, while when z ^ zn the diffraction affects the field 
shape. 



3 Lorentz- Gauss beams 

So far we have evaluated the paraxial propagation of the beam essentially by 
using an integral approach. In fact, once we knew the field on the source plane 
Vq{x, y) we were able to obtain it on every plane as dictated by Fresnel theory, 
namely we had 

V{x,y,z) = --^^^^ I jj,{Ul)eM{^Yz^^-S.f + {y~vfmd^ (25) 



5 



Actually we used (|25|l in Fourier space to obtain the plane wave spectrum on 
z and from it, by mean an inverse Fourier transformation, we pointed out the 
field V(x, y, z). If we let 

V{x,y,z)^cx^[ikz)](x,y,z) (26) 

where /(a;, y, z) is the slowly varying part of y, 2:), we know that, if T^(a;, y, z) 
fulfils the integral H25I) . then f{x,y,z) fulfils a differential equation, known as 
paraxial wave equation, that in tridimensional space holds 

V|/(x,y,z) + 2zfc/,(x,y,z) = (27) 

where is the two-dimensional Laplace operator in the transverse plane, i.e. 
^t/ = fxx + fyy and fz is the partial derivative respect to z variable. This 
equation has very special properties, and in particular we here are interested 
to its symmetry features. In fact there exists a beautiful theory, due to the 
mathematician Sophus Lie, that allows to perform an analysis on the symmetry 
groups associated to a particular differential equation (or, more in general, to 
a system of differential equations). Essentially, the theory says that there exist 
some differential transformations which act like operators on system's solutions 
and that lead to others solutions of the same equation, when applied to an 
existing and known starting solution. The theory bases itself upon a certain 
number of theorems in the context of differential geometry, and we will not 
enter in details; however, in Appendix we put a proof of the derivation of the 
symmetry group we utilize in the following (see, for example, in llOj for the 
theory foundations). This theory was used, for example, by Wunche JT] to 
show that Hermite- Gauss and Laguerre-Gauss beams (with complex argument) 
can be generated from the fundamental Gaussian beam simply by applying to 
it the powers of certain differential Lie operators. Among all the symmetry 
tranformations associated to l|27|l we concentrate upon only one which states 
that if f{x, y, z) is a solution of the aforesaid equation also will be the function 
f*{x,y,z) where 

r*f X 1 r _f!_±i!_i f( ^ y 

J [x,y,z) ^^^^/^exp[ ^2(1^^^/^)1/4^ -^z^' 1^-^/^' 

(28) 

with wq and L two real parameter. The property in 1)28(1 was recently proved, 
utilizing two different approaches, by Kiselev that has utilized a separation 
variable method and by Gutierrez- Vega et al. [S]), in which the authors obtained 
the same result by inserting a well-constructed trial function into paraxial wave 
equation. It is interesting to note that the way we indicated here is a third 
different method, a method which possesses the merit to underline which is the 
foundation of this result, i.e. an hidden symmetry and which shows that all 
beams with a Gaussian envelope are intimately connected to the paraxial wave 
equation. As a consequence of (|28|l we immediately conclude that also exist the 
Lorentz-Gauss beams, namely 

TT*/ X ^TJ"^ exp(ifcz) . , , X z . x z 

^^y^^'l + iz/L' l + iz/L^ ^ ^l + tz/L' 1 + iz/L^^ 
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The shape of this kind of field depends by the choices of the parameters Wx,Wy, wo 
and L. In particular if we choice wq < Wx,Wy, we obtain a beams that behaves 
like a Lorentz beam near to the z axis and like a Gaussian beam far from it. If 
we put L ~ zji, Lorentz-Gauss beams for wq — w^iwy) shows a field shape 
more defined around z axis with respect to the Lorentz beams (I18|l . i.e. on the 
X — y plane the LB has a wider bell-shaped field function with respect to (|29f) . 

4 Conclusion 

A new kind of tridimensional, rectangularly-symmetric, realizable scalar optical 
beams has been introduced. On the source plane these beams are the product 
of two indipendent Lorentz function and the exact analytical expression for 
the field on a generical plane, under paraxial regime, has been derived. In 
particular it is interesting to note that it deals of a rare case in which one is in 
the presence of physically realizable fields, the propagation of which is known 
in closed-form, that does not possess a Gaussian envelope term. This kind of 
field can represent a valid candidate to modelize the shape of fields generated by 
certain laser sources, as double hetero junction (DH) Gai-^AlxAs lasers. Using 
a Lie group approach we introduced also the solution obtained by modulating 
the aforesaid beams with a Gaussian envelope. 
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A Symmetry groups of scalar paraxial wave equa- 
tion 



Consider the 3-D scalar paraxial wave equation 

\7'^u{x,y,z) + 2ikuz{x,y,z) = (30) 

We wish to show that this equation admits a class of solutions that are mod- 
ulated by a two-dimensional Gaussian envelope. First of all, to this equation 
it is possible to associate a manifold M E X x U, where X = 5Rp, with p = 3 
meaning number of independent variables and C/ = 3?'', with q ~ I meaning the 
number of dependent variables. On such a manifold it is also possible to define 
a tangent vector field in the following form 

V = ^{x, y, z, u)dx + r]{x, y, z, u)dy + r(x, y, z, u)dz + (j){x, y, z, (31) 

In the context of Lie group theory tangent vector fields are the generators 
of symmetry transformations through the following relationship 

V'(e,v) =exp(ev) (32) 

where if) is the transformation induced by the vector v and e is a real parameter 
characterizing the group. A symmetry transformation is a map that allows 
to pass from a starting point {x,y,z,u) on the manifold M to another point 
(x' ,y' , z' ,u') on the same manifold by mean the relation 

{x', y', z', u') = ?A(e, v)(x, y, z, u) = exp(ev)(x, y, z, u) (33) 

To find the explicit expression of the vector field in JSU, one can utilize the 
following procedure. It is necessary to extend the space X x U \n order that it 
also contains the second order derivatives, an operation said prolongation. By 
consequence the prolongated tangent vector field becomes 

V = id^ + iidy + Td, + c^du + + c^ydu^ + c^'du^ + + 

+ ct^yydu^^ + r^du^^ + r^du^^ + r^du^^ + (34) 

where we have dropped the dependence from the variables x,y,z,u. We can 
rewrite it as 

p 

i=l J 

by defining J — (ji, j2, ■■■ji), 1 < jz < P j 1 < I < n where n is the equation 
order. All the coefficients in (|34|l are expressible in terms of 6 ''' ^-iid (p and 
their derivatives through the formvila 

p p 

<j,' ^Dji<p-J2cuj,r)+J2euj,^ (36) 

i=l 4=1 

where Ui — du/dxi and uj^i — duj/dxi with Xi generical variable and Dj 
representing the total derivative. Among all such 13 coefficients only <^^^, (p^y, (f)^ 
are important for our purpose. Indeed, under certain hypotheses, it is possible 
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to obtain all the symmetry group of transformation of scalar paraxial wave 
equation through the condition 

^[uxx + Uyy + 2iku^] = (37) 

that, by taking into account relation (|34|1 implies 

+ + 2ik4>' = (38) 

On performing the calculations as in 13611 we find 

0"^"^ = (jixx + Ux(2(j)xx - S.xx) + Uy{-rixx) + U^{-Txx) + ul{(j)'''^ -2£,uu) + 

~\~ UxUy(^ 2tIxu) ~^ '^x'^zi '^'^xu) ^" ^xx'^xi 3^-;^) + Uxx'^yi Vu) ~t" 

+ UxxUz{-Tu) + UxyUxi-2T]u) + Uxx{(f>u - 2^a;) + wiJ(-3^„„) + 

+ Uyx{-2rix) + Uzxi-'^Tx) + UxU^xi-'^Tu) + ul.Uy{-r]uu) + 

+ u1u^{-Tuu) (39) 

f/)™ = (l)vy + Uy{2cj)yu - r^yy) + Ux{-£,yy) + U:^{-Tyy) + U'^y{cj)^'' - 2'qyu) + 

H~ 'Ux'^yi^ '^^yuj H" 'Uy'^zi^ '^'^yu) ~^ '^yy'^y(~^'^u^ H~ '^yy'^xi^ ^"u) 

+ UyyU^{-Tu) + UxyUy{~2S^u) + Uyy{(f>u " 2r]y) + u^y{-iriuu) + 

+ Uyx{~2^y) + Uzy{~2Tu) + UyU^y[-2Tu) + ulUx{-^Uu) + 

+ ulu^{-Tuu) (40) 

= (l>z + Uz{2(j)u - T^) + Ux{-£,z) + U:^Ux{-^u) + Uy{-ri^) + 

+ UzUy{-riu) + ul{-Tu) (41) 

Now we equate the right and left-hand side homologous terms appearing in H38|l 
and finally we obtain 



Oyy = <Pt 



(42) 

2(02;« - ixx) - Cyy = "Ct (43) 
'^{^yu ^xx) ^xx ^y (44) 

(45) 
(46) 
(47) 
(48) 

0n« - (49) 

where we have defined an auxiliary variable t = z/{2ik). Solving this system 
of equations is not difhcult but we report here the result only, that can be 
verified by substitution, 

r = ci + 2c2t + Acst^ (50) 

— C2X + Ac^tx — 2C4< + C5 (51) 

fl = C2y + ^c-ity - 2c4i + cg (52) 



2Cx = 


Tt 


2riy = 


n 


'^x '^y 


Tu = 




= 
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(j) = [-€3(2;^ + y^) + C4(.x + y) - 4c3t + cj]u + a{x, y, t) (53) 

where a is a generical function and Ci{i = 1,2, ...7) are integration constants. 
Each of these constants is related to a particular generator of symmetry 
which can be obtained by letting tidily all the constants to zero except the ith. 
Among all such a transformations there is the following one 

V = itxd^ + Uydy + U'^dt + {-x^ - - it)udu (54) 

From H33() - (|54|l follows that, if we indicate as u{x, y, s) a solution then also the 
following one represents a valid one 

u*(x,y,s) — - ; — exp[— e- — — r^-rl^f-; : — , -; — —, — , -; : — ) (55) 

^ ^ l + 4es (l + 4es)^ H + 4es 1 + 4es' 1 + 4es^ ^ ' 

On coming back to the old coordinate z and letting e/2k — 1/L and Wq = 1/e, 
it holds 



1 x"^ + y"^ 

1 + iz/L '^^'^^^wlil + iz/L)'""-l + iz/V l + iz/L' l + iz/L' 

(56) 

that represents the result we were looking for. 
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Figure 1: Modulus of field's amplitude of a Lorentz beam on: a)sourcc plane, 
i.e. z = h) when z = O.lzu c) when z = O.^zu d) when z = zr. The 
parameters are choosen as follows: A = 0.6328/im, Wx = Wy = w = lO^A. 
Zr = 1-nw^l\ = 3.976m 
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